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Abstract. We present a simple model for the evolutionary states of neutron stars in pre-low-mass X-ray binaries 
during the main sequence stage of donor stars. It is shown that for typical parameters some of the neutron stars 
in these systems can not accrete matter from winds of their companions in a stable way. Accretors are found if 
neutron stars have magnetic fields about 10 13 G and higher and/or in close systems with orbital periods about a 
few days and shorter. 
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1. Introduction 

The increased ability to observe very faint X-ray 
sources with modern X-ray observatories (XMM-Newton, 
Chandra) prompted new interest in systems with low ac- 
cretion rates ~ 10 10 -10 12 g s _1 like isolated neutron stars 
(NSs) or NSs in detached binaries. Among isolated NSs 
the fraction of accreting objects is not expected to be 
high, and accretors are expected to be dim (see Popov 
et al. 2003 for a review and discussion) . 



Recently IWillems fc Kolbl l|2003|) discussed evolution 
and observational appearance of pre-low-mass X-ray bi- 
naries (pre-LMXBs) (see also Belczynski & Taam 2003). 
Pre-LMXBs by definition represent an evolutionary state 
prior to LMXB. These systems are detached binaries in 
which NSs accrete matter from the stellar wind of its com- 
panion and have detectable accretion luminosity. 

Willems and Kolb found two maxima in the distri- 
bution of accretion luminosity of NSs in pre-LMXBs at 
~ 10 31 and ~ 10 28 erg s _1 . The second maximum cor- 
responds to very small accretion rates similar to those of 
isolated NSs accreting the interstellar medium. One can 
suspect that for such small mass flow a NS can avoid ac- 
cretion at all. IWillems fc K olb ( 2003|) did not explore the 
evolutionary state of NSs in pre-LMXBs assuming, that all 
of them can accrete. Here in a very simple model we try to 
calculate if NSs in pre-LMXBs can reach the stage of ac- 
cretion, or they are on the propeller or ejector stages (see 
Lipunov 1992 for detailed description of different stages 
and Popov et al. 2003 for a short introduction). 



2. The Model 

We assume that the evolution of a NS starts at the stage of 
ejection, then continues to the propeller stage (including 
subsonic propeller, see Davies & Pringle 1981) and Anally 
a NS can reach the stage of accretion. Stages are sepa- 
rated by critical periods (see detailed derivation of them 
in Lipunov 1992). A NS is born as an ejector and starts 
to slow down. The ejector stage ends when P = Pe: 



Pe = 2tt 



2^ 2 



C 4 FMacc 

cP/2tt 



1/4 



10/4£V 1 / 4 V 6 1/2 s,JZi <Rg. 



Here R\ 
R G = 2GM NS /V 2 



4.8 10 9 P cm - light cylinder radius, 
= 3.7 10 14 Vg~ 2 cm - radius of accre- 
tion, = fj,/10 30 G cm 3 - magnetic moment of a NS, 

V 6 = FTotal/lO^mS" 1 , Vrotal = V^wind + ^orb) 2 , " " 

number density of the stellar wind at the NS position. In 
all formulae here and below we assume NS's moment of 
inertia to be equal to 10 45 g cm 2 . 

Note that the condition R\ < Rq is not necessary ful- 
filled in the binaries we are going to study. We discuss this 
topic in the Sec. 4. 

Subsonic and supersonic propellers are separated by 
the condition of equality of the co-rotation radius, R co = 
(GMP 2 /47T 2 ) 1 / 3 , and radius of the magnetosphere, Ra- It 
corresponds to a critical period: 



P ss = 2 5 / 14 n(GM ] 



NS 



.-5/7 



(^/Macc) 3/7 



« 300 ^ V 3/ V 6 9/7 s, R A < R G . 

For some systems the relation Ra < Rg can be violated. 
We discuss it below in the Sec. 4. 
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Accretion starts when P — Pbr (unless a NS is not on 
a georotator stage). 



Table 1. Parameters for the models in fig. 1 



F br = 610V3o /21 



M" 5/7 m- 4 / 21 s. 



(1) 



Here M s — M ajcc /10 s gs~ 1 - accretion rate, m = 
Mns/Mq. We use numerical coefficient 6 10 6 in eq.Q fol- 
lowing Davies & Pringle (bu t scaling it to t he value of 
M acc typical for pre-LMXBs). Ilkhsanovl l|2003|) suggested, 
that it should be a factor 7.5 larger. If it is correct, then 
the time scale of the subsonic propeller stage is longer. 
However, this stage is not well understood, and we try to 
be conservative, and whenever possible choose values to 
lower timescales of pre-accretion stages. 

The output of our calculations consists of two 
timescales: ejector time scale, t&, and propeller timescale, 
tp (see Prokhorov et al. 2002, Lipunov 1992 for more de- 
tails). 

At the ejector stage a NS is spinning down ac- 
cording to magneto-dipole formula (approximatelly P = 
3 10 _4 ^3oiyrs )• Propeller stage consists of two distinct 
phases: supersonic propeller and subsonic propeller (see 
Davies & Pringle 1981 and Ikhsanov 2003), so we define 



^supcr "T" 



Actually, there should be also an inter- 



mediate stage between supersonic and subsonic propeller, 
when turbulence becomes supersonic above the magneto- 
spheric radius, but inside Rq- Spin-down during this 
intermediate phase should be similar to the subsonic pro- 
peller rate, so they are united here. Supersonic propeller 
is a relatively short stage (normally shorter than each of 
ejector and propeller stages). For supersonic propeller we 
choose a very efficient spin-down (see Popov et al. 2000). 

Taking all together we can write the following equa- 
tions for durations of the ejector stage (iE), supersonic 
(isuper) and subsonic (i su b) propeller stages: 



fE = 0.810V 3 o 1 '^ 1/2 ^yrs, 



t. 



-8/7 



n 



super — 1.3 10 /I3Q 

t S ub = 10 3 /%, 2 m P b r yrs 



3/7T/9/7 



yrs, 



(2) 

(3) 
(4) 



Eq. i s taken fol lowing IPopov et all ll200Cl). ea. <I31 - 
fromlshakural 11197.51) . and eq. (0J - from lDavies fc Pringle! 
lll98lh . 

To calculate orbital velocity, V rb, we assume circular 
orbits and equal masses Mns = M* = 1.4 Mq (Willems 
and Kolb discuss Mq < M* < 2 Mq). In our calculations 
the orbital velocity is always smaller than the stellar wind 
velocity 

To calculate the accretion rate we use the following 
equation (see Willems & Kolb 2003): 



-A/acc. — 
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To 



2 2 
r_ 



1 + q R* 

2/3 2 a 



-3/2 



M*. 



Here P* - stellar radius, q = Mns/M*, M* - the rate of 
stellar wind mass loss. We used several different values for 





M acc , 


M30 


Torb , 
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Mod. 


M/r, yr~ 


O 




days 


km s — 




num. 


1(T 13 


3.8 10 10 


0.1 


1 


660 


0.03 


1 


icr 13 


3.8 10 10 


1 


1 


660 


0.04 


2 


10- 13 


3.8 10 10 


10 


1 


660 


0.04 


3 


10- 13 


5.2 10 8 


0.1 


30 


660 


0.008 


4 


10- 13 


5.2 10 s 


1 


30 


660 


0.01 


5 


10- 13 


5.2 10 s 


10 


30 


660 


0.01 


6 


10- 13 


5.0 10 6 


0.1 


700 


660 


0.003 


7 


1 n-13 
10 


5.0 10 6 


1 


1000 


660 


0.003 


8 


nr 13 


5.0 10 6 


10 


1000 


660 


0.003 
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Table 2 


. Parameters for the models in fig. 
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M„ 


M acc , 






Vwind , 


t Hupor 

'sub 


Mod. 




1 gs" 1 




days 


km s _1 




num. 


10" 13 


3.3 10 11 


1 


1 


330 


0.07 


1 


10" 13 


3.3 10 11 


10 


1 


330 


0.08 


2 


10- 13 


7.6 10 9 


1 


30 


330 


0.02 


3 


10- 13 


7.6 10 9 


10 


30 


330 


0.03 


4 


HT 12 


3.3 10 12 


1 


1 


330 


0.13 


5 


10- 12 


3.3 10 12 


0.1 


1 


330 


0.1 


6 


10- 12 


7.6 10 10 


1 


30 


330 


0.04 


7 


lO" 12 


7.6 10 10 


10 


30 


330 


0.05 


8 


lO" 12 


2.9 10 10 


1 


1 


1320 


0.03 
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M* and K, ind = (3 ^/2GM*/R* (see tables, 330 km s" 1 
corresponds to = 0.5, 660 km s _1 - to (3 — 1 and 1320 
km s _1 - to j3 = 2). Note, that actual accretion rate can 
be lower than this value due to heating, magnetospheric 
and hydrodynamical effects (see discussion in Perna et al. 
2003 and Popov et al. 2003). 

To derive many equations in this paper we will use the 
following relation which parametrises the accretion rate 
by n and V^otai: 



M 



= 710 11 

Af»/(47ra 2 K, ind ) 



n V 6 3 gs" 



density of matter, m P 



n = pm p , p 
- proton mass. 

Fo r M* on the main sequence stage IWillems fc Koltj 
l)2003h used the value 10" 13 M yr" 1 . We follow them in 
that choice and also check higher value 10~ 12 Mq yr , 
which was also discussed by this authors. Definitely for 
smaller M* all timescales (£e, ^super, ^sub) become longer, 
for higher M* - shorter. 

Mass loss can be increased significantly when a donor 
star leaves the main sequence, but we do not follow these 
stages in our analysis. Also we do not discuss the possibil- 
ity of a formation of a binary system due to capturing of 
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a companion (in that case the picture can be completely 
different, as far as an old evolved NS can be captured). 

We want to note once again that our model is a very 
simple one as far as we want just to illustrate the impor- 
tance of the effects of magnetorotational evolution of NSs 
in pre-LMXBs. Discussion on possible complifications can 
be found below in the Sec. 4. 

3. Results 

Our main results are presented in the figures. Each point 
represents a model with particular set of parameters (see 
tables). The three lines correspond to the total time of 
evolution prior to the stage of accretion (ixotai = + tp) 
equal to 1, 5 and 10 billion years. 

Note, that tp oc ig, a w 1.3. This happens because 
of the following reason. In eq. P|l tw. oc /i _1 n -1 / 2 and 
VTotaJ is n °t varied significantly in our modeling as far as 
Kvind > V OI b- If we substitute eq. JIJ into eq. 10} , then 
we have i su b oc /i~ 26 / 21 n~ 5 / 7 . (Remember, that in tp we 
have i supor << t su ^ [see tables].) As we see we have nearly 

tp OC <g 3 . 

Vertical line in the figures corresponds to tp, = 
2.3 10 9 yrs - the main sequence lifetime of a star with 
M* = 1.4 Mq. (Roughly, the main sequence lifetime 
can be estimated as: logiMS = 9.9 — 3.8 log(M*/MQ) + 
log 2 (M^/Mq), see for example Lipunov 1992.) If for some 
reason the propeller stage appears to be very short (for 
example, subsonic propeller is not operating) then NSs to 
the left of this line can accrete. 

As one can see most of models from table 1 have 
^Totai > 2.3 10 9 yrs. It means that these NSs can not ac- 
crete. Accretion is more probable onto strongly magne- 
tized NSs as far as their spin-down timescales are shorter 
on both ejector and propeller stages. To allow accretion 
with a rate ~ 10 10 g s" 1 a NS with B ~ 10 12 G should 
be slowed down to very long periods ~ 10 5 s which takes 
a long time. 

On the second figure we present models with txotai < 
5 10 9 yrs. These are systems with strongly magnetized NSs 
or/and with shorter orbital periods or/and with smaller 
wind velocity. 

4. Discussion and conclusions 

Our calculations w ere made for constant magnetic fields . 



As it was shown bvlColpi et al.1 1 1998h . iLivio et alJ v ,. 

IPopov fc Prokhorov ( 200fJ) realistic parameters of decay 
in the case of isolated NSs normally make the propeller 
stage longer. Parameters of the magneto-rotational evo- 
lution of isolated NSs and NSs in pre-LMXBs are very 
similar (for example, there is no huge accretion to speed- 
up the process of field decay), so the results for the former 
ones can be applied to the latter (and vice versa). In that 
sense our estimates of tp and tp for B — const should 
be lower limits if the field is decaying. For decayed fields 
most of NSs in pre-LMXBs can be expected to stay on the 
propeller stage. For high Kvmd a NS instead of accretion 



V) 



M 
O 




Log t E , yrs 

Fig. 1. Crosses with numbers represent different set of 
parameters of NSs (see table 1). Three curves corresponds 
to frrotal = ts+tp equal to 1 Gyr (dashed), 5 Gyrs (dotted) 
and 10 Gyrs (solid). The vertical line corresponds to tp = 
2.3 10 9 yrs. 
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Fig. 2. Same as in fig. 1 for parameters in table 2. 



can stay on the so-called georotator stage (Lipunov 1992). 

It happens for Vxotai ^ 470 /%) n. 1 / 10 kms -1 (see Popov 
et al. 2003). 

In this short illustrative note we neglect several effects, 
which can change characteristic periods and radii or/and 
spin-down rates. Most of the time scales used above can 
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be considered just as rough estimates (see discussion in 
Popov et al. 2003). However, we believe that qualitatively 
our results are valid as far as our main aim is only to 
illustrate importance of the effects of magneto-rotational 
evolution of NSs in pre-LMXBs. 

Let us briefly discuss possible violations of the rela- 
tions i?i < P G , i?A < P G (for details see Lipunov 1992). 
We start with the relation between R\ and Rq. For high 
magnetic fields Pe becomes large. It necesserely means 
large Ri. On the other hand K,i n d is also large, which 
means small Rq. For systems 3, 6, 8, 9 in the table 1 and 
for system 9 in the table 2 R\ > Pg at the moment of tran- 
sition from ejector to propeller. In that case equations for 
Pe and £e are the following: 

P E2 «180^ 3 n- 1 /3 1 /- 1 /3 S] (5) 
iE2«3.510 11 M3o 4/ V 1 /3^/ 3 yrs . (6) 

For large velocities and magnetic fields normally > 
£e2> so the stage of ejection is shorter for Pi > Rq (com- 
pare eqs.[21inj). For our investigation it is crucial to know 
if this effect is strong enough to let the system to cross the 
line t~E — 2.3 10 9 yrs on the figures. Only for the model 9 
from the table I it is so (in the figures all data point are 
plotted for Pi < Pg)- For all others this shortening of £e 
is not strong enough to change the conclusions made in 
this short note. Also we have to mention that if a system 
leaves the ejector stage after the violation Pi < Pg then 
the accretion rate is expected to be very low. 

Now let us discuss the relation between Pa and Pg- 
Often Pa < Pg (for example it is so for all observed 
bright X-ray binaries, for isolated NSs which can become 
accretors etc.). For systems with highly magnetized NSs 
and large wind velocities it is possible, that Pa > Pg- It 
is important for example to estimate the critical period 
P ss and t su b, isupcr- Alfven radius is determined by the 
following equations: 

„ ( ^ y /7 _,. 

V2M acc% /GlW 
8.210V3oV 2/7 lf /7cm > Pa < Pg, 

V M acc V* J 

7.610 u A4o 3 n~ 1/6 ^" 1/3 cm» R A2 > Ra- 
in the later case the equation for the critical period is: 



we do not take this effect into account. Study of georo- 
tators in binaries is a separate (and probably promissing) 
subject. At that stage the re is a possibility of a particular 
"magnetic accretion", see iRutledgel l|2003|) . 

In the discussed type of binaries there are possibilities 
for episodes of accretion. As far as low-mass stars can have 
episodes of eruptive activity it is possible that NSs can 
accrete matter due to fluctuations in the stellar wind. Post 
main sequence stages are relatively short, so no significant 
spin-down is expected. There is just a quick increase of the 
accretion rate: increase of M» and decrease of Vmnd both 
work to increase M a cc- Non-accreting systems, which are 
situated in t^-tp-plane close to the line — 2.3 10 9 yrs 
(like model 1,5,9 in the fig. 1) can start to accrete if M acc 
is increased when the star leaves the main sequence but 
the system is still detached. Systems with highly eccentric 
orbits can have episodes of accretion in periastron. 

We conclude, that as in the case of isolated NSs many 
NSs in pre-LMXBs (at least on the main sequence) are not 
at the stage of stable wind accretion. Accretion is possible 
when the donor is on the stage of main sequence if a NS 
is strongly magnetized or in close binaries. 
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P ss2 = 2n V4 - 2.5 V 1/2 s. 

For systems with Pa > Pg durations of the propeller 
substages are different from eqs.(J21EJ. These systems later 
appear not as normal accretors but as georotators and so 



